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problem of twist elasticity of fluctuating ribbons using the simplest worm 
like chain model (WLC) remains elusive. In this paper, we approach this 
problem by taking first a mechanical approach and then incorporating ther- 
mal effects in a quadratic approximation applying the Gelfand-Yaglom (GY) 
method for computing fluctuation determinants. Our study interpolates be- 
tween mechanics and statistical mechanics in a controlled way and shows 
how profoundly thermal fluctuations affect the elasticity of semiflexible poly- 
mers. The new results contained here are: 1) a detailed study of the minimum 
energy configurations with explicit expressions for their energy and writhe 
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fluctuations around the local minima of energy and approximate analytical 
formulae for the free energy of stretched twisted polymers derived by the 
Gelfand Yaglom method. We use insights derived from our mechanical ap- 
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of thermal fluctuations. From the derived formulae, predictions of the WLC 
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merical simulations and experiments. 

Keywords Twist Elasticity • Gelfand Yaglom • Thermal Fluctuations • 
Worm like chain 

Supurna Sinha 

Raman Research Institute, Bangalore, India 
Tel.: +91-80-23610122 
Fax: +91-80-2361 0492 
E-mail: supurna@rri.res.in 

Joseph Samuel 

Raman Research Institute, Bangalore, India 



2 



PACS 64.70.qd • 82.37.Rs • 45.20.da • 82.39.Pj • 45.20.Jj • 87.14.gk 
1 Introduction 

In recent years, the theoretical study of the elasticity of semifiexible poly- 
mers has emerged as an active area of research. These studies are motivated 
by micromanipulation experiments 3,6, 14,7 on biopolymers in which sin- 
gle molecules are stretched and twisted to measure elastic properties. These 
experiments are designed to understand the role of semifiexible polymer elas- 
ticity in the packaging of these polymers in a cell nucleus. Twist elasticity 
plays an important role in several biological functions. DNA, a long molecule 
(between microns and meters in length) which carries the genetic code, is 
neatly packed into the tiny cell nucleus, just a few microns across. The first 
step in packaging DNA in a cell nucleus involves DNA-histone association 
which makes use of supercoiling in an essential way. The process of DNA 
transcription can generate and be regulated by supercoiling|43] . 

The single molecule experiments have been theoretically analyzed using 
the popular Worm Like Chain model, which gives an excellent account of the 
bending elastic properties of DNA 24 . When the twist degree of freedom 
is explored by rotationally constraining the molecule the theoretical analysis 
gets complicated by subtle topological and differential geometric issues. In a 
previous work 37 we have treated these subtleties in the language of differ- 
ential geometry and topology. Our present study builds on earlier mechanics 
based research in this area 12,4,13,7 and goes on to incorporate the effects 
of thermal fluctuations in a quadratic approximation. We evaluate the parti- 
tion function and present an explicit analytical form which can be compared 
against measured elastic properties or numerical simulations (29) . 

Two parallel streams of research have evolved in the study of semifiexible 
polymer elasticity. Some researchers use classical elasticity which is concep- 
tually simple, but ignores thermal fluctuations which are important to the 
problem. Others use a statistical mechanical approach which properly takes 
into account thermal fluctuations and entropic effects. This latter approach 
is beset by the geometrical and topological subtleties we alluded to earlier. 
In the real biological context of a cellular environment, these polymers get 
constantly jiggled around by thermal fluctuations and therefore a statisti- 
cal mechanical approach is necessary. Biological processes involve the entire 
range of flexibilities: Actin filaments, microtubules and short DNA strands 
are energy dominated, while long DNA strands have appreciable entropy. The 
quantitative measure of flexibility is L/Lp, where L is the contour length of 
the polymer and Lp its persistence length (16/i for Actin and 50nm for 
DNA). In the energy dominated short polymer regime classical elasticity 
gives a fair understanding of the problem. In the absence of external forces 
and torques, long polymers require a statistical mechanical treatment since 
the entropic contribution is not small. In this paper we consider polymers 
subject to a stretching force and torque and develop a theory of thermal fluc- 
tuations around the minimum energy configurations. These minimum energy 
configurations come in two families, the straight line family (SLF) and the 
writhing family (WF). We present a treatment based on classical elasticity 
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and then incorporate thermal fluctuations around the minimum energy con- 
figurations. Members of the writhing family are far from straight and so our 
analysis goes far beyond perturbation theory about the straight line. While 
our general analysis applies to polymers of all lengths, some of our analytical 
formulae are derived in the long polymer limit. Our treatment is analytical 
and complementary to numerical studies [26, 45, 291. 

Two extreme situations are relatively well understood. If the polymer is 
nearly straight, one can use perturbation theory [2"71H01rT5] about the straight 
line to calculate its elastic properties. In the opposite extreme when the poly- 
mer is wrung hard, it buckles and forms plectonemic structures P51[TB1[T5] 
which are stabilized by the finite thickness of the DNA. The transitional 
regime where the polymer is neither straight nor plectonemic is not as easy. 
Perturbation theory about the straight line is not applicable and neither do 
we benefit from simplifications that arise from the energy dominated plec- 
tonemic regime. More specifically, let the stretching force F be 2ksT/Lp or 
higher. As we turn the bead to twist or "wind up" the molecule, it initially 
remains approximately straight. At high links the polymer buckles and winds 
around itself to form plectonemic structures. The intermediate link regime, 
in which the polymer is neither straight nor plectonemic is the subject of 
this paper. (These three regimes of low link, intermediate link and high link 
have also been referred to in [55] as straight line, buckling and supercoiled 
regimes.) We will study the straight line and the writhing solutions and fluc- 
tuations around these solutions. 

The central goal of this paper is to derive explicit analytic expressions for 
the free energy and the writhe of a twisted stretched polymer. Derivatives of 
the free energy expressions computed here connect to measurable quantities 
which can be probed via single molecule experiments or simulations. We 
arrive at the free energies by focusing on the stationary points of the energy 
functional dominating the partition function. While there is interest in the 
elastic properties of DNA at all length scales, the experiments [3,43, 14 deal 
with long polymers, many times the persistence length of about 50nm. Many 
theoretical treatments [30, 2 , 33 , 29 also focus on this regime of long polymers, 
which leads to some simplification. 

In order to formulate the problem we start from an idealized experiment. 
Let us imagine that a twist storing semiflexible polymer (like dsDNA or 
filamentary Actin) is attached to a glass slide at one end and a bead at the 
other, so that its tangent vector at both ends is constrained to point in the z 
direction. Supposing the glass slide end to be at the origin, we apply a force 
F and torque r (both in the z direction) on the bead. The elastic response 
of the polymer can be described by its extension z = r.F (the position of 
the bead) and its link (the number of times the bead has turned about the z 
axis). Real experiments on DNA differ from the one described above mainly 
in that they deal with long polymers and it is usually the Link rather than 
the torque that is held constant. Though the experiments [14,21 could easily 
be adapted by using a feedback loop to maintain constant torque. 

Actually, we will be interested in a wider class of thought experiments 
of which the above is a particular realization. We may, for instance, use 
more general boundary conditions, fixing the tangent vector at the ends to 
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t(0) = ti and t{L) — if. (The force is still along z and the torque along 
if.) As theorists, we may also explore the entire range of flexibilities and 
parameters, which are not practical in the currently accessible regime. Our 
study is not restricted to long polymers, though we will sometimes specialize 
to this case for simplicity. We will also regard all our parameters (contour 
length, persistence length, force, torque, temperature, elastic constants and 
boundary conditions) as tunable without worrying about how this may be 
experimentally achieved. In the mechanical limit our analysis reduces to the 
classical theory of beams and cables. Our larger aim is a comprehensive un- 
derstanding of the entire range of parameters covering the range from beams 
and cables to DNA. Thus we are not immediately concerned with any one 
experiment, but a class of them. Detailed treatments of actual experiments 
are contained in Rcfs. [25,7 . For example these treatments take into account 
the chiral nature of the DNA molecule, which we ignore. We are more con- 
cerned here with mathematically working out the predictions of the simplest 
worm like chain model rather than a detailed modelling of specific experi- 
mental data. However, our focus is experimental relevance and we calculate 
from the model, quantities like link-extension relations which are measured 
in laboratories. 

The paper is organized as follows. Sec II deals with the mechanics of 
semiflexible polymers. We start by noting that the problem of computing 
the partition function of a twist storing polymer separates into two parts, a 
simple Gaussian over the twist and a harder problem involving the writhe. 
We derive the Euler-Lagrange (E-L) equations describing the minimum en- 
ergy configurations. The solutions of the E-L equations can be expressed 
using elliptic integrals. Sec III deals with the role of thermal fluctuations. 
We perform the second variation of the energy and explain the methods we 
use to compute the thermal correction to the free energy due to fluctuations 
about the classical solutions. This gives the main result of our paper, an 
approximate expression for the free energy of a stretched twisted polymer 
which takes into account thermal fluctuations around nonperturbative solu- 
tions. We then suggest calculational schemes (Sec. IV ) for a more ambitious 
treatment of thermal fluctuations. We end with some concluding remarks in 
Sec. V. 



2 Mechanics 

We model the twist storing polymer by a ribbon (x(s), e J (s)), (i = 1,2,3) 
which is a framecQ space curve. x(s) describes the curve, i(s) — its 
tangent vector and e l (s) the framing, s is the arc length parameter along the 
curve ranging from to L the contour length of the curve. x(0) = since 
one end is fixed at the origin. The tangent vectors at both ends t(0),t(L) are 
fixed. A force F = Fz along the z direction is applied at the free end at x(L) 
and also a torque r along the fixed end tangent vector t(L). We suppose 



1 Note that we do not use Serret-Frenet framing, which is popular in this field. 
Serret-Frenet framing becomes ill defined over locally straight pieces of the curve. 
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83 = t(s) and refer to e±(s) as the ribbon vector that describes the twisting 
of the polymer about its backbone. Defining the "angular velocity" 17, via 

#± = « x e, (1) 
as 

The "body- fixed angular velocity" components of H are J?j = e$ ■ fi. The 
expression for the energy of a configuration is 

£ (C)= [ -[A(f2l + nl)+Cf2l)ds- { F.tds - 27rrLk, (2) 

where Lk is the number of times the bead is turned around the if axis. The 
mathematical problem we face is to compute the partition function 

Z (F,t)=5>p-!^. (3) 
c B 

In Eq. ((3]) , the sum is over all allowed configurations of the ribbons [37] . The 
ribbon can be closed with a fixed reference ribbon that goes a long way in 
the if direction, makes a wide circuit and returns to the origin along the ii 
direction^. In the calculations below, we will set both A and fc^T equal to 
unity and restore them when necessary. 

We use the celebrated relation [46ll5] decomposing the link into twist and 
writhe: 

Lk = Tw + Wcw (4) 

where Tw — 2t: J Qs(s)ds and Wcw is the writhe. The writhe is a non-local 
quantity defined on closed simple curves: Let the arc length parameter s 
range over the entire length Lq of the closed ribbon (real ribbon + reference 
ribbon) and let us consider the curve x(s) to be a periodic function of s with 
period L . Let R(s,ct) = x(s + a) — x(s). The Calugareanu- White writhe is 
given by[MlMIMTl] 

1 f L " f L °- ,dR(s,a) dR(s,a), - 
Wcw= A ds d<r[^iK^i].£ (5) 



Air J J Q+ ds da 

Because of White's theorem Eq.@, we find that the problem neatly splits 
16,27,37 into two parts 

Z q (F,t) = Z 1 (F,t)Z 2 (t) (6) 

where 

Zi (F, r)= J C[x(«)] exp -fx [x(s)] (7) 

and 

f L 1 di di f L 
fi[x(«)] = / ds - — ■ 5- - / F • tds — 2ittWcw- (8) 
J 2 ds ds J 

2 The reference ribbon is supposed to be nowhere self intersecting or south 
pointing |37| . So we also suppose 7^ —£. 
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Z%{t) is given by a Boltzmann sum over framings with £2 given by 

£ 2 (t) = J ds[—Ql - t2ttTw] (9) 
Z2 (t) is easily evaluated as a Gaussian integral and gives 

Zi{r) = exp [^] = exp [-G 2 {r)\. (10) 

Thus, the problem that remains is to compute Z\(F,t) (Eq.0 ) which de- 
pends only on the curve x(s) and not its framing. 

This problem is hard because of the appearance of the writhe Wcw[x(s)], 
which is a non local function of the curve. However, we will make progress 
by noting that variations of the writhe are local[37|l5B]. We will compute the 
partition function assuming that for high stretch forces, the sum over curves is 
dominated by configurations near minima of the energy. The approximation 
consists of using an expansion of the energy about the minimum energy 
configuration and keeping fluctuation terms about the minimum to quadratic 
order. This section deals with the minimum of energy and the following one 
with fluctuations. We thus derive explicit analytical expressions for the Gibbs 
free energy of the polymer Q (F, r) which can be used to compute its elastic 
response. 

Theoretically, it is easiest to deal with the constant torque ensemble. In 
the quadratic approximation, the conjugate ensembles are equivalent [2011421 
I551I4T] even for polymers of finite length. This equivalence of ensembles holds 
exactly for long polymers, since these are at the thermodynamic limit. (See 
however, the remarks below concerning stability in different ensembles.) 

In this section, we ignore thermal fluctuations and consider a purely clas- 
sical mechanical analysis to study the configurations of a torsionally con- 
strained stretched semiflexible biopolymer like DNA. While such treatments 
exist in the literature 16,31,22,4,12,13 , we will present a slightly different 
perspective based on analogies with classical mechanics, which help us to 
incorporate thermal fluctuations. 

The central quantity of interest is the bending, stretching and writhing 
energy, ^ ^ 

£{C) = \ [ lids - [ F.tds - 2ittW C w (11) 
^ Jo Jo 

of a space curve x(s) whose tangent vector is i(s) = The tangent vector 
is varied subject to the boundary conditions t(0) = U, i(L) — if fixing the 
tangent vector at both ends of the curve. Remembering that the variations 
of writhe are local [371136] . we arrive at the Euler-Lagrange equations (see 
Eq.dMJ) below) from Eq.(ITT)l: 

- 1 - F = r(i x i) - 7*, (12) 

where the term jt arises since t ■ t = 1. The problem is formally similar 
to a symmetric top, a fact that was well known to Kirchoff The analogy is 
useful for integrating the Euler Lagrange equations. We use quotes for the 
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analogous top quantities. The "Kinetic energy" is given by T = ht-t and the 
"potential energy" is V = F • i. The total "energy" 

H = T + V = -i ■ i + F • t = ^(9 2 + sin 2 9tp 2 ) + F cos 9 (13) 
is a constant of the motion as is the z component of the "angular momentum" 

J z = (t x i) z — ri z = sin 2 9ip — t cos 9, (14) 

where we have introduced the usual polar coordinates on the sphere of tan- 
gent vectors. Using these "constants of the motion" we reduce the problem 
to quadratures as described in [TJ5]. The basic equations are 

§2 v j? a (Jz+t cosO) 2 

— = H — Fcos9 5 15 

2 2 sin 2 6 V ' 

{Jz+TCOsff) 

</>= — a ■ (16) 

sin 

Setting u = cos 9, we find that 

u = ef(u) (17) 



where f(u) — y/V(u),e — ±1 and 

V(u) := 2{H - Fu)(l - u 2 ) - (J z + tu) 2 (18) 

is a cubic polynomial in u. Integration of Ea. (|17[) gives u(s) in terms of elliptic 
integrals and Ea. (IT5|) gives <p(s). Further integrating the tangent vector gives 
us the polymer configuration x(s) in real space. Since our polymer is under 
stretch (unlike that of Ref. [11] which is under compression) we assume F 
is positive. V(u) is positive for large positive it, and negative (or zero) at 
u = 1 and u — — 1. In order for there to be a physical solution there must 
be two real roots, b, c (turning points) in the physical range of u — cos 9: 
— 1<c<6<1 and a third real root a > 1. We have ordered the three roots 
so that a > b > c. The boundary data must lie between c and b. c > Ui,Uf < b 
. The "motion" goes from m to Uf, possibly passing through turning points 
on its way. 

The integration constants H and J z are determined by the length of the 
polymer and the boundary conditions. From Hamilton- Jacobi theory we can 
write the energy of the solution of length L going from Ui,tpi to Uf,iff as 

£(ui,Uf,<pi,tpf,L)= duf(u) -HL + J z $ (19) 

J Ui 

where L and = iff — (fi, given by 

f Uf du , , 

L= 77^ 20 

Jui f( u ) 
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and 

•-L m htu) (21) 

determine H. and J z and /i(u) = (J 2 +ru)/(l — it 2 ). These integrations must 
include turning points if any. For instance, if u goes from ttj to Uf via c, 
= XT* /c f ' a ll° wm g f° r sig n changes in the integrand Ea. (fP7|) . 
We will use these relations below in computing the fluctuation determi- 
nant in the next section. The expressions for <£, L and £ are elliptic integrals, 
but as these functions are not as widely known today as they were a cen- 
tury ago, it is more useful and appropriate to leave them as integrals. It is 
quite easy to numerically evaluate these integrals on a computer and plot 
any function which may be of interest. 



2.1 Symmetric Boundary Conditions: 

As an illustrative example, we specialize to boundary conditions ij = if = z. 
These boundary conditions imply that 9 = is a point on the solution. Since 
9 has to be finite and non- negative at 9 — 0, we have J z = — r and T~L > F. 
The form of V{u) simplifies to 

V{u) = (1 - u)[2(H - Fu)(l +u) - r 2 (l - u)}, 

and Eq. p^|) reduces to 

cp= — . (22) 

We want to find the minimum energy configurations subject to the writhe 
constraint. These configurations satisfy the Euler-Lagrange equations. The 
simplest solution is the straight line t(s) = z or 9(s) = for all s. This 
configuration globally minimizes £{C) 

£ ST = -FL (23) 

but since Wcw = 0, it cannot accommodate writhe. For this we need the 
"writhing solution" , in which the tangent vector t starts from z at s = and 
9 increases to a maximum of 9q at s = so, the turning point and then returns 
to z at s — 2sq. The period of the orbit is P = 2sq. As was shown in [36] , 
the configuration is a local minimum of the energy only if the length L of the 
polymer is equal to P (and not a multiple of it). 

It has been shown earlier (36) that the local minima of the energy are 
"good curves" [30] i-6 they satisfy the conditions of Fuller's [T^IITS] theorem 
and so the writhe can be computed from a simple local formula Q 

Wcw _ r- *,(! ,)»(,) 

Below we will sometimes drop the subscript on Wcw, it is understood that 
we only deal with "good curves" . 

3 As a side r emark we mention that our approach is guided by geometric phase 
ideas 39,23,34,35 
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2.2 Long Polymers 

One gets further simplification by confining to the limit of long polymers. In 
this limit, the "period" P (= L — > oo). As a result the doubly periodic func- 
tions that appear for finite L become simple trigonometric functions (since 
one of the periods tends to infinity) . We will briefly restrict our attention to 
infinitely long polymers in this section to arrive at an explicit analytical form 
for the "writhing family" of solitons. To reach this limit, the length Eq. (|2TH) 
must diverge. For % > F, the function V{u) = f 2 (u) vanishes linearly as 
u —¥ 1. The corresponding integral converges and does not lead to L —> oo. 
Thus we must require that V(u) vanishes quadratically as u — > 1 so that 
Iu JJu) diverges logarithmically at the upper limit. This situation obtains 
if V(u) has coincident roots at the upper limit. Now, demanding that n^r~n 

vanishes as we take u 1 in Eq. (TT71) gives us the condition T-L = F. So the 
form of V{u) simplifies to 

V(u) = (1-u) 2 [2F(1 + u)-t 2 } 

The turning points of u are at u max — 1 and u m i n — uq, where 

The solutions are found by elementary integration [21[r2] 

u(s) = (1 — uq) tanh 2 n(s — sq) + uq (25) 

ip(s) — ^— — — — Arctan[— ta,nh/i(s — sq)], (26) 

2 r 

where /i 2 = F — r 2 /4. The energy of the writhing family parametrized by 
F, r and L is given by 

£ w = -FL + 8^ tanh ^— - 2tttW C w{t). (27) 

Since we have already established 36 that the writhing family are "good 
curves" |30) we may compute the writhe using the simpler "Fuller formula" 
Wf, which (apart from normalization) measures the solid angle swept out 
by the unique shorter geodesic connecting the tangent vector to the north 
pole. We find by a straightforward calculation that the writhe of the writhing 
family is 

2signumr .2u , ixL. . , 

W C w = W F = arctan[-4tanh^-], (28) 

7T \t\ 2 

in agreement with [7]. Including the twist energy Eq. (fTU|) — t 2 L/(2C), the 
"Gibbs" free energy Q (F, r) of the writhing family in the mechanics approx- 
imation is 

G${F,t) = -FL + 8Aitanh[^] - 2t:tW C w{F,t) - Lt 2 /{2C) (29) 
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Fig. 1 Gibbs free energies Q{F,t) (from a purely mechanical approach) of the 
writhing family (thick red curve) and of the straight line family (thin blue line) 
plotted against torque r for L = 10, F = 2 and twist elastic constant C — 1.4. 
A, ksT, Lp are set to 1. Note that the writhing family has higher energy and that 
this difference narrows as one nears buckling. 

where Wcw is given by Ea. (|28|) . Correspondingly, the "Gibbs" free energy 
Q (F, t) of the straight line family in the mechanics approximation is 

gf L (F,r)=-FL-^ (30) 

These two Gibbs free energies are plotted in Figure 1 as a function of r 
for a fixed force. In our expressions Eg. ([281291) for the energy and writhe 
of the writhing family, we have retained expressions like tanh[fj,L/2], which 
become 1 in the long polymer limit provided /i is not small. If one drops these 
expressions our energy difference coincides with [12J after correcting for an 
overall factor ( y/F/2) , a presumed misprint in [12] . 

The writhing family always has a higher energy |12| than the straight line 
family for fixed force and torque. However, this difference is not extensive 
in the length L and remains finite as L — > oo. The energy difference at zero 
torque is 8V FA. At higher torques (see Fig. 1) this energy difference narrows 
and disappears at the buckling torque (v AFA, for long polymers). At the 
buckling torque, the writhing solution merges with the straight line and there 
are no classical solutions beyond this torque. 

The writhing family becomes relevant when there is a need to accom- 
modate writhe. As one can see by toying with a tube, applying Link to the 
ends of the tube results in deformations which cause the tube to deviate 
considerably from the straight line. If one applies a fixed link to the tube, 
the elastic response of the tube depends on C, the twist elastic constant. In 
Fig. 2,3, we plot the relative extension £ = z/L versus link for a fixed force 
for large and small values of C. If C is large (Fig. 2), there is a critical link 
above which the SLF ceases to exist. The WF is then the only solution to the 
mechanics problem. Further increasing the link causes the writhing family to 
buckle into plectonemes. Thus the writhing family represents the transition 
from the straight line to plectonemes. 
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Fig. 2 Relative extension £ of the polymer versus applied link for the straight 
line family (thin blue line) and the writhing family (thick red line) for a stretching 
force of F = 2 (in dimensionless units of ksT /Lp) and L — 10. This plot has a 
twist elasticity C = 10 in dimensionless units. At high Link the straight line family 
ceases to exist; only the writhing family can store writhe. 




Fig. 3 Relative extension £ of the polymer versus applied link for the straight line 
family (thin blue line) and the writhing family (thick red line). This plot uses C=1.4 
and a stretching force F = 2 and L = 10. Similar curves are seen in simulations 

mi- 

For smaller values of C, the behavior is as shown in Fig. 3, which plots the 
relative extension of the polymer versus link. If one takes a straight ribbon 
and twists it, it remains straight and acquires link, following the thin blue 
line in the figure. At a critical torque, the straight line meets the writhing 
family (thick red line in the figure) and becomes unstable, follows the red 
line of the writhing family and buckles. This is exactly the behavior seen 
in the simulations of Ref.[2Jj]. See the section of the curve from bo to b\ in 
Figure 2 of [H] . The rest of the curve in this figure of Ref . [35] is not relevant 
to our analysis since it involves self contact of the polymer and goes beyond 
our analytic approach. A drop in extension at buckling has also been seen in 
the experiments of [14) . 
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3 Fluctuations 

In a thermal environment, the polymer fluctuates around its minimum energy 
configuration. This contributes a fluctuation term to the partition function 
and gives a thermal correction to the classical energy computed in the last 
section. The elastic properties of the polymer are profoundly affected by 
these fluctuations. Interpretation of the experiments on DNA elasticity [33] 
requires an inclusion of thermal effects. By taking appropriate derivatives of 
the free energy computed here one can find the theoretical predictions of the 
worm like chain model for the experimentally accessible extension versus link 
or torque twist relations. (See |14j for experimental graphs of these.) 

For a system with a finite number of degrees of freedom, q l , i — 1..N with 
a potential V(q), the mechanical energy in equilibrium is given by the value 
of the potential V(q*) at its minimum q* where 



where we formally write the gradient of V as dV and the Hessian matrix as 
ddV . (The Hessian is positive at local minima. ) By performing a Gaussian 
integral, we find that the free energy in the presence of thermal fluctuations 
is 



Our objective here is to present a calculation of the free energy of twisted 
scmiflcxible polymers by considering fluctuations about the mechanical so- 
lutions of the last section. At fixed force and torque we have the minimum 
energy classical solution, which satisfies 5£b s — 2tttSWcw , where Wcw is 
the writhe and £bs is the bending and stretching energy. Formally, what we 
wish to do is to expand the energy and writhe around the classical configura- 
tion retaining terms to quadratic order. The quadratic integral is a functional 
Gaussian integral which can then be performed and gives an answer in terms 
of the determinant of the fluctuation operator. The problem reduces to com- 
puting the determinant of the fluctuation operator. 

Computing the functional determinant of the WLC polymer involves some 
subtleties. If one evaluates the determinant as an infinite product of eigenval- 
ues, one finds that the product does not converge and needs regularization. 
This is an artifact of the fact that we have modeled the polymer as a space 
curve (or ribbon), which has an infinite number of degrees of freedom. Phys- 
ically, we know that the polymer consists of a finite number of atoms and 
has finite free energy. There is a well developed mathematical theory to deal 
with such situations. Multiplicative (additive) constants can be dropped in 
evaluating the partition function (free energy) and the infinities do not affect 
physical predictions. Another subtlety arises from zero modes of the classical 
solution. If the energy functional (and boundary conditions) has continuous 
symmetries which are broken by the classical solution, there are "Goldstone 
modes" which cost no energy and have zero eigenvalue. Such zero modes need 
separate handling. We will face both these problems below. 



dV\ q . -0 
ddV\ r > 



(31) 
(32) 



T = Vict) + -k B T\ogdetddV\q*. 



(33) 
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Let us now consider a family of polymer configurations parametrized by 
a. The first variation of the energy is (where ' denotes the a derivative). 

£'[*(*)] = / -i'-F- 1') ds (34) 



W'[*»] = 4 / (i(s)xi(s))-i'(s)ds. (35) 



and the first variation of the writhe is 

f 

The Euler-Lagrange equations read £' — 2ittW' = or 



[ {-i-F-T(i(s)xi{s))}-i'ds = (36) 
Jo 



for arbitrary t' satisfying t' ■ i = 0. 

Now consider the second variation of the energy 



£" [*(«)] = / i-t - I F-i"+ I t -t ds (37) 



1 f L A , . 1 ' L 



and the writhe 

W"[f(s)] = — / (tx i) ■ t" ds + — I (t(s) xt )- i' ds. (38) 
2tt Jo 2-7T J 

Computing £" — 2tttW" we find that the t" terms combine, thanks to the 
Euler - Lagrange equations Eq. flT!?]) . to give 7 Li' ■ i' ds. (Use the identity 

i-i" +i' -f = o). 

The form of the second variation functional is 

£"-2tttW"=[ {t ■ t - T(t cl (s) x t) ■ t' + j e i(s)t' ■ i'} ds (39) 



This is a quadratic form in t' which we assume is non negative (since we 
suppose that t c i (s) is a minimum of the energy, either local or global) . The 
subscript on t c i(s) and 7 c z( s ) indicates that it is the classical solution about 
which we wish to compute the determinant of this quadratic form. Varying 
t' to find the fluctuation operator (analogous to ddV above) we find the 
eigenvalue equation (for i' satisfying t c i.t'(s) = 0.) 

6i' = -i - Tt ci x i' + ~{ cl i' = \i' (40) 

which gives us the spectrum of the fluctuation operator. If we could find all 
the eigenvalues of O, we could compute its determinant. While this works in 
special cases (like the SLF) it not usually possible in general. 

To compute the determinant of the fluctuation operator O, we use a 
technique due to Gelfand and Yaglom (GY) [TT1II25] . We consider the energy 
Eg. (1111) as a function of the boundary conditions at the end of the polymer. 
(This is similar to considering the Action in classical mechanics as a function 
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of the end points in Hamilton Jacobi theory.) We consider classical solutions 
which go from t\ at s = to tf at s — L. Although our physical problem 
fixes ti and tf, the calculational technique (GY) requires that we consider 
variations of the energy under such changes. The result of GY is that the 
computation of an infinite dimensional functional determinant (the determi- 
nant of the operator O) reduces to the computation of a finite dimensional 
one. The final answer is expressed in terms of the variation of the classical 
energy with respect to variations of the boundary conditions. Let £{ti,tf,L) 
be the energy of the configuration of length L that goes from ti to tf (where ti 
and t f are the small variations of the boundary conditions about the actual 
boundary conditions in our problem). The infinite dimensional functional 
determinant after regularization is reduced to the computation of a simple 
finite dimensional determinant: 

detO = [det-l^r 1 (41) 
dtidtf 

evaluated at the physical boundary data. From Hamilton- Jacobi theory we 
can write the quantity in square brackets in Ea. (|4TI) as 

dpf 



det-f (42) 



where is the "final momentum" p(L) of the classical configuration of 
length L that reaches tf from U. 

The method of GY is used in many fields of physics (see for example [5]) 
for computing fluctuation determinants. Ref.|llj contains an application of 
the technique to polymer physics in which it is used to compute the thermal 
corrections to Euler buckling of polymers under compression. Our present 
application is to stretched and twisted polymers. The great advantage of 
the GY technique is that one does not have to diagonalize the fluctuation 
operator O to compute its determinant. All one has to do is to linearise the 
Euler Lagrange equations and study the Jacobi fields. In many cases, it is 
not practical to diagonalize the fluctuation operator, while computing the 
Jacobi fields is a much easier task. 

We saw in the last section that the minima of the energy come in two 
families, the straight line family and the writhing family. Since the GY tech- 
nique may not be familiar to many readers, we illustrate its use by applying 
it to a case where the answer is known, the straight line family. We calculate 
the fluctuation determinant by two methods, first explicitly diagonalizing the 
fluctuation operator and then by the GY technique. We then apply the GY 
idea to the writhing family to obtain new results. 



3.1 Straight line Family 



For the straight line family, i c ;(s) = z (for any F, r). As is well known [l2l 
[27] . there is a range of (F, r) (r < r c — \JAF + (f;) 2 )) for which the straight 
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line family is stable against perturbations. The eigenvalue equation Eg. ([40)) 
reads i 

- i - t(z x t ) + Ft' = At', (43) 

where t'.z — The eigenvalues are easily worked out and the thermal correc- 
tion to the free energy computed. This involves regularizing a divergent sum 

^b^E^i logi^ 2 + I 7^ _ )' which can be done by standard methods. The 
final answer turns out to be detO = (sinh/xL//x) 2 which leads to the thermal 
correction to the free energy 

, „, r sinhLH n 

k B T log [ (44) 

A* 

in agreement with Ref.|18j. 

This computation can also be performed by line arising the Euler La- 
grange equations about the straight line. After some simple transformations, 
it reduces to the inverted oscillator (in a magnetic field) and the solutions can 
be written in terms of hyperbolic functions (unlike the trigonometric func- 
tions of the usual oscillator). Varying the relations between final and initial 
data, we easily find 

5Z(L) = (exp-irL/2)[(cosh M L)5Z(0) + At - 1 (sinh^L) ( 5P(0)] (45) 
6P(L)= (exp-iTL/2)[(i(smhiJ,L)5Z(0) + (coshiJ,L)5P(0)] (46) 

where Z = t x + it y are complex coordinates in the tangent plane at the 
north pole and P the conjugate momenta. Setting 5Z(L) — and solving for 
SP(L), performing the differentiation Eq. (|4"2l yields the determinant detO = 
[det §§mj ] = (sinh/iL/^i) 2 as before. Computing k B T j 2 log det\ §^jy| 2 gives 
us the same answer as Eq. (|4"4"|) . 



3.2 Writhing Family 

For the writhing family, the eigenvalue equation Eq. (|40[) is considerably more 
complicated (note the presence of the functions t c i(s) and "f c i{s)) and the 
eigenvalues and eigenvectors are not available in closed form, but we can still 
use the GY method. The writhing family breaks the azimuthal symmetry 
(rotation about the z axis) that is present in the energy and the boundary 
conditions. As a result the solution has zero modes which need careful han- 
dling. There are two routes open to us and we briefly describe both of them. 
We can alter the boundary conditions so that they break the azimuthal in- 
variance. For instance, we could make at least one of ti and tf not point along 
the z direction. The other route is to exploit the symmetry of the rotational 
invariance and treat the azimuthal degrees of freedom exactly. This leaves a 
one dimensional problem (9 or it), which we can treat either approximately or 
exactly, depending on the desired degree of accuracy. Since the calculations 
are involved if straightforward, we outline the method and present the final 
answer. 
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3.3 Asymmetric Boundary Conditions 

We suppose that Ui,(fi and Uf,ipf, the initial and final tangent vectors are 
not both in the z direction. We make a variation in the initial positions, while 
holding the final positions fixed. We are interested in the resultant variation 
in the final momenta. To compute this, we proceed as follows exploiting the 
integrability of the system. We will keep the final configuration (uf,ipt) fixed 
and perform variations in (ui,(pi) the initial configuration and look at the 
variations in <P, L and £ . These are 

° L = 77 — \ + "557 5 ™ + "57 ° J 
f{ui) d% oJ 

S<P = -j-^r h( Ui ) + — 5H + — SJ 
f(ui) dH dJ 

58 = €if(ui)5ui + 5J$ + J5$ - 5UL - HSL. (47) 

The Jacobian of interest is Eq. (|4"2"]) . which can be expressed as (writing pf 
for p UJ and J z for p Vf ) 

d(p f ,J z ) = d( Pf ,J z ) d(H,J z ) 
d(ui,cpi) d(H,J z ) d(ui,<pi)' 



From 



we find 



e /% =^ /)/( l-4) = 7^)-^^ 



Ksjj-y 7 ^ n r [si = A [5J Z ) ^ 




where A is a 2 x 2 matrix defined by Eq. 
Setting SL = and 5<P = 5(tp 2 - tfi) = —S(fi in Ea.([4"?|) we find 

w-v-/- c (^ (50) 

where the 2x2 matrices V and C are defined in Eq. ([5"0"f . From 

(detO)- 1 = det {^A , (51) 

we finally arrive at 

detO = (det^)- 1 (detP)(dctC)- 1 

= /(«,)/(«,) det 2? (52) 

rdLd$ d<PdL, 

-^M[dndJ^mdj ]e ^ (53) 
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where V is the 2x2 matrix defined in Eq. (l50l) . The final answer for the total 
free energy including thermal corrections is 

Q{F,r,k B T) = E cl + *f log [«*,/(«,)/(«,) [f||§ - §|§]] (54) 

These derivatives are evaluated at constant it,, it/. While this final answer 
may appear to be abstract, it is in fact quite amenable to numerical methods. 
All the functions appearing here can be expressed Eq. (|19l2QI21j) in terms of 
elliptic integrals and numerically evaluated. 

Although Ea. (|54p is written for asymmetric boundary conditions, the 
thermal correction to the free energy for long polymers with symmetric 
boundary conditions can be teased out of it, since the free energy for long 
polymers is expected to be independent of the boundary conditions. We take 
the limit of Ui — > 1, which leads to /(it,-) vanishing. In Eg . ([54]) the terms that 
multiply f(ui), we need only keep those terms that diverge as itj — > 1. This 
results in a simple expression for the leading correction to the free energy: 

AQ W = k B TLn = k B TL^F - t 2 /4 (55) 

In Ea. (|55p we have dropped lower order terms which are not extensive in 
L. Notice that the numerical value of the thermal correction to the free 
energy is the same as for the straight line family. This can be understood 
physically since for any finite /x, most of the polymer is straight. We would 
not in general expect the fluctuation energies of SLF and WF to be equal in 
short polymers. The expressions Eqs. (|53M55I) contain the main results of this 
study; an analytic calculation of the thermal corrections to the free energy. 

Our main result expressing the fluctuation determinant in simple terms 
sheds considerable light on the stability of polymer configurations. Note that 
the determinant vanishes as either itj or Uf approaches a turning point. This 
signals the appearance of null eigenvalues for the Hessian or Fluctuation 
operator Eq. (|4T)|) . Such effects have been well studied in optics and mechan- 
ics pQ, where they reveal the appearance of caustics, focusing and conjugate 
points. If the polymer configuration goes from u, to it/ passing through turn- 
ing points, it is not stable against small perturbations. The Gibbs free energy 
G(t, F) can be lowered at constant torque by such perturbations. At constant 
torque, such configurations will decay to lower energy configurations. 

However, if one works in the corresponding Helmholtz ensemble, at con- 
stant link, which is how most experiments are done, the allowed perturba- 
tions are also required to maintain constant writhe (supposing for simplicity, 
C to be infinite) and this instability disappears. However, if m and it/ are 
separated by two turning points, there are two independent perturbations, 
which lower the Energy. Working to second order in the perturbation, one 
of these lowers the writhe and the other increases it. One can take a linear 
combination of these two perturbations to maintain constant writhe. Such 
configurations are unstable at both constant torque and constant Link. Note 
that polymer configurations that extend for more than a period P arc there- 
fore unstable to perturbations in accord with the findings of |36j . They are 
not local minima of either the Helmholtz or Gibbs energies. 
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4 Towards Statistical Mechanics 

Till this point we have relied on the mechanical approach, dealing with in- 
dividual polymer configurations and considering small fluctuations around 
them. This has the advantage that we know which individual configuration 
we are dealing with and can show [3 6) that the dominant configurations admit 
a local writhe formula. We now move towards a more statistical mechanical 
approach and integrate over configurations. In the process, we lose touch with 
individual configurations, but we gain by getting a fuller treatment. From the 
mechanical approach, we have gained confidence that the sum over config- 
urations is dominated by "good curves" that admit a local writhe formula. 
Using this idea, we will now extrapolate to more general situations where 
the thermal fluctuations are appreciable and continue to use a local writhe 
formula. There is a logical jump involved here and we have argued elsewhere 
[57] that this leads to a good approximation provided that the polymer is 
under stretch. 



4.1 Symmetric Boundary Conditions 

If we assume that the boundary conditions are azimuthally symmetric, we 
can exploit the symmetry to treat the tp degrees of freedom exactly by "in- 
tegrating them out". We reduce the problem to a single variable 9 (or u). In 
the process of reduction the effective potential picks up an extra contribu- 
tion from the thermal fluctuations of the tp degrees of freedom. As a result 
the effective potential explicitly depends on the temperature. The origin of 
this effect is that path integration in general coordinates has subtleties com- 
ing from the measure^, as pointed out by Edwards and Gulyaev [10]. These 
techniques have been applied [32] to compute the quantum corrections to 
soliton energies. Our calculation is mathematically similar, but physically 
different since we deal with thermal corrections rather than quantum ones. 
Our treatment below closely follows these earlier studies, with appropriate 
modifications. 

After integrating out the ip degrees of freedom, the 9 motion is governed 
by an effective potential 



k B T ,(1- cos 6>) k B T 

Veff{6) = Fcos(9- -2— +r 2 -i -f (56) 

m ; SL P 1 + cos 9 8L P sm 2 9 



and the earlier equation Eq. (|17[) is replaced by 

W 2 = V T (u) = (/tW) 2 = 2{{n+ k ^)-Fu)(l-u 2 )-r\l-uf+ k ^- (57) 

Note the explicit temperature dependence in Ea. (|5fJ]) . Note that u = 1 is 
no longer an allowed solution. The straight line is destabilized by entropic 
effects. Indeed it would be hard to distinguish between SLF and WF here 

4 This can also be seen by operator methods, where it emerges from the need to 
make the reduced differential operator self adjoint. 
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Fig. 4 Relative extension £ versus link for a fluctuating polymer for the straight 
line family (thin blue line) and the writhing family (thick red line). The inner curves 
correspond to F — 2 and the outer ones to F = 15. 



because we are no longer dealing with single configurations but ensembles of 
them. However, we can still continue to use a local writhe formula for torques 
large enough that the turning point c is positive. Such curves have tangent 
vectors entirely in the northern hemisphere and are "good" curves. 

We now treat the 9 (or u) motion approximately by considering the "con- 
figurations" u(s) which satisfy the Euler-Lagrange equations. One finds the 
"classical" contribution by evaluating the energy on the "classical" solution. 
("Classical" is in quotes because it is a bit of a misnomer, since the effec- 
tive potential includes contributions from thermal fluctuations.) The thermal 
correction to the free energy can now be computed: 

dL 

F(T, F, t) = £„ cv , (T, F, t) + k B T/2 log [ - 4(7* - F) ( — )] (58) 

In the first term in Eq. (|58[) £»ci" is computed using the thermally corrected 
Vt{u)- In the second term which is already of order fcgT, we have used the 
old V(u) since the difference is higher order in ksT. L is given in terms of % 
by formula Ea. ((20|) . These formulae are numerically tractable (analytically 
they are still elliptic integrals) and can be used to work out the predictions 
of the model. 



4.2 Statistical Approach 

In fact, it is possible to also treat the u degree of freedom exactly but with 
more work. The effective potential Eg. (155)) gives rise to a Schrddinger equa- 
tion 

-l^ + Veff{e)i> = ^. (59) 

Solving the Schrddinger equation gives the full propagator, which in turn 
gives us the free energy of the polymer. It is convenient to pass to imaginary 



20 



"time" which changes the relative sign between the kinetic and potential 
terms. This turns the potential over V e ff — > —V e ff. For r 2 > the po- 
tential has an attractive well at 9 — tt coming from the last two terms in the 
effective potential (Eq.lSB)). The potential near 9 — tt diverges as — f /{tt— 9) 2 . 
As a result, writhe can be stored at arbitarily small energy cost. This pathol- 
ogy of the model is caused by paths that wind round the south pole and has 
to be regulated by imposing a cutoff pi[T5] . This is a reflection in this model 
of the plectonemic transition 0. As these authors point out, a cutoff corre- 
sponding to the effective diameter of the DNA molecule is necessary. In fact 
the twist experiments reported in [5] show the polymer making transitions 
between the plectonemic phase and the normal phase. It would be interest- 
ing to compare the predictions of this theoretical model with the experiment. 
Possible points of contact are the relative time spent in each phase and the 
mean rate of transitions, both of which can be calculated from the model. Of 
course a comparison with the real experiment will also involve the dynamics 
of the bead as done in [7J. Our study deals purely with the polymer. 

For long polymers, we can use the eigenstate expansion with states ip n 
and energies e n 

K(t h tf, L) = E n ip n (ii))ip n (if) exp [ - e„] (60) 

is dominated by the ground state and compute the free energy as the ground 
state energy eg of the Schrbdinger operator in Eq. (|59[) . 



5 Conclusion 

We have given a detailed analysis of the mechanics and fluctuations of writhing 
polymers under stretch. Our main results are 1) an explicit characterization 
of the mechanics of the writhing family, its energy and writhe 2) a study 
of fluctuations about these equilibrium configurations and formulae which 
can be used to calculate the thermal correction to the free energy. In par- 
ticular, we have presented explicit approximate analytic expressions for the 
free energy fEas. (|M] ) and of a semiflexible polymer ribbon. We do not 
have room here for a detailed elucidation of the "phase diagram" which re- 
sults from the free energy expressions above. However, our work shows that 
approximate analytical treatments of the problem of writhing DNA are possi- 
ble. Such analytical treatments are complementary to simulations (numerical 
experiments) and real experiments. 

A curious feature of the energy (and free energy) of the writhing family 
is the cusp at r = 0. On differentiation this results in a discontinuity in the 
link between positive and negative values of r. The discontinuity is visible in 
Figure 2 (the gap in the thick red curve) and is also present in simulations 
(see Figure 2 of [1H]). This feature can be traced to the fact that the r = 
member of the writhing family has self intersections and there is a consequent 
jump in the writhe across it. 

We started this paper with a mechanical approach and moved progres- 
sively to a more statistical mechanical point of view. The advantage of me- 
chanics is that we are dealing with single configurations of the polymer and 
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can therefore check whether the local writhe formula applies. The work of [36 
shows that it does apply to all the local minima of energy (at least for a class 
of solutions (J z = — t), which includes symmetric initial data). We expect 
from here that small fluctuations about these "good curves" are also "good 
curves" . Further extrapolating this idea, we expect that the present models 
based on a local writhe formula are a good approximation to the exact, but 
difficult to solve models based on the non local writhe. This result supports 
the idea proposed in |37] that in a statistical sense, one can approximate self 
avoiding models by south avoiding ones. 

We mentioned before that in mechanics, the straight line is incapable of 
accommodating writhe. This is shown in Figs 2 and 3, where the extension of 
the thin blue line is maximal independent of the applied link. This situation 
changes drastically when thermal fluctuations are incorporated as in Fig. 
4. The straight line family no longer has maximal extension (the thin blue 
line curves down). Thus the applied link is stored entirely in the thermal 
fluctuations of the molecule. The molecule has helical fluctuations and the 
molecule stores writhe by preferentially writhing with one helicity. This is in 
sharp contrast to the writhing solutions, which can store writhe even without 
thermal fluctuations. 



A mechanical treatment is expected to work well for short polymers. 
For long polymers without applied force, the system is entropy dominated 
and therefore our treatment does not apply. For long polymers under a very 
high stretching force (for example the regime F « 10 in L TJ corresponding 
to lpN), the entropic fluctuations are tamed and a mechanical approach is 
again possible. As the force is lowered, entropic effects begin to gain impor- 
tance. Our treatment considers fluctuations about non perturbative solutions 
of the writhing family and therefore goes beyond perturbation theory about 
the straight line. We have argued elsewhere [37] that the free energy predicted 
by "south" avoiding models are a good approximation to that of "self avoiding 
models" provided the polymer is under stretch. This is certainly true at low 
link, since perturbation theory about the straight line applies. At extremely 
high links, the same is true for a different reason: the configurations in which 
the polymer winds around itself in the self avoiding model are mimicked by 
south winding configurations in the south avoiding model. At intermediate 
links, we expect that the polymer is in a coexistence phase |25j with some 
of its length storing the writhe in the plectonemic phase and the rest of it 
at a torque below the buckling torque. Our study suggests that even in this 
regime which is far from perturbative, the two models give approximately 
similar predictions. As mentioned earlier, this is based on our demonstration 
[36] that the dominant configurations (the local minima of the energy) are 
good curves. Only a simulation will convincingly demonstrate how low one 
can go in stretching force before the approximation breaks down. We hope 
to compare these theoretical predictions with computer simulations in the 
future. 
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